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Abstract: The mechanism of the chemical reaction of H 2 0 with three stabilized Criegee 
intermediates (stabCI-OO, stabCI-CH3-00 and stabChc-OO) produced via the limonene 
ozonolysis reaction has been investigated using ab initio and DFT (Density Functional 
Theory) methods. It has been shown that the formation of the hydrogen-bonded complexes 
is followed by two different reaction pathways, leading to the formation of either OH 
radicals via water-catalyzed H migration or of a-hydroxy hydroperoxide. Both pathways 
were found to be essential sources of atmospheric OH radical and H2O2 making a 
significant contribution to the formation of secondary aerosols in the Earth's atmosphere. 
The activation energies at the CCSD(T)/6-31G(d) + CF level of theory were found to be in 
the range of 14.70-21.98 kcal mol -1 . The formation of a-hydroxy hydroperoxide for the 
reaction of stabCbc-00 and H2O with the activation energy of 14.70 kcal mol 1 is 
identified as the most favorable pathway. 

Keywords: volatile organic chemicals (VOCs); ozone; ab initio methods; ozonolysis 
reaction mechanisms; limonene 
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1. Introduction 

The stabilized Criegee intermediates (stabilized carbonyl oxides) from the ozonolysis of alkenes can 
react with various atmospheric compounds [1-5], particularly with the formaldehyde, H 2 0, NO x , SO2, 
H2SO4 and CO and many others. One of the most important reactions in the atmospheric chemistry is 
the reaction of stabilized Criegee intermediates with water, the dominant constituent of condensable 
vapors in the Earth's atmosphere. This reaction is considered to be one of the major degradation 
reactions in the Earth's atmosphere [6,7]. It is also known that this reaction leads to the formation of 
a-hydroxy hydroperoxides (HOCH 2 OOH or HMHP), organic acids, ketones, aldehydes, OH radicals, 
and H2O2 [6,8-19]. These species have been detected in both the ambient air and precipitation, in both 
forested and urban areas under polluted conditions [15,20-25]. HMHP may act as an enzymatic 
inhibitor of peroxidases, while H2O2 is a very important oxidant [26-28]. The H2O2 contributes to acid 
precipitation via the conversion of SO2 into H2SO4 [29], which is the key atmospheric nucleation 
precursor. The reaction of stabilized Criegee intermediate with water can be an additional essential 
source of OH radicals [30]. The oxidation of H 2 0 by stabilized Criegee intermediates in ozonolysis of 
alkenes can also contribute to the formation of secondary organic aerosol (SOA) [31,32]. 
Jonsson et al. [31] corroborated the aforementioned conclusion for the ozonolysis of limonene, 
A -carene, and a-pinene at low concentrations of the aforementioned species and noticed that the 
particle number and mass concentration increases with increasing relative humidity The reaction of 
stabilized Criegee intermediate with water is also important for the water purification and waste water 
processes, where ozonolysis is widely used [33]. 

Previous experimental and theoretical studies on the reaction of H 2 0 with CH 2 00 and CH3HCOO, 
(CH3)2COO and the stabilized Criegee intermediate from ozonolysis of limonene indicate that this 
reaction is the three-step mechanism [2,6,34-38]. At the first step, a hydrogen-bond complex is formed (a). 
At the second step, two decomposition pathways for the hydrogen-bond complex are possible: 
water-catalyzed hydrogen (H) migration would lead to the formation of OH radical via 00 bond 
breaking (b). At the final stage, the a-hydroxy hydroperoxide is formed by the addition of the water 
molecule to the stabilized Criegee intermediate. Sauer et al. [14] found that a-hydroxy hydroperoxides 
RR f C(OH)OOH formed in the reaction of the stabilized Criegee intermediate and water may have 
different fates depending on the chemical nature of the R and R f substitutes (see Scheme 1). 

Being produced by more than 300 species, limonene or 4-isopropenyl-l-methyl-cyclohexene, is the 
most abundant monoterpene in the Earth's troposphere [39,40]. The reaction of ozone with limonene, 
which has both endocyclic and exocyclic double bonds, is one of the most important oxidation 
processes in the Earth's troposphere. Ozonolysis of limonene is an essential source of OH radicals, 
H2O2 and plays a significant role in the formation of atmospheric aerosols via the oxidation of H2O by 
stabilized Criegee intermediate [31]. There are four stabilized Criegee intermediate formed from 
cleavage of the 03-Limonene primary ozonide: stabCI-00 (See stabCI-00 structure in Figure 1) and 
stabCI-CH3-00 (See stabCI-CH3-00 structure in Figure 2) formed from the endocyclic primary 
ozonide decomposition and stabCbc-00 (See stabCbc-OO structure in Figure 3) and stabCH 2 00 
formed from the exocyclic primary ozonide decomposition [41]. 
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Scheme 1. Different fates of RR f C (OH)OOH formed in the reaction of the stabilized 
Criegee intermediate and water depending on the chemical nature of the R and 
R f substitutes. 
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The mechanistic diagrams for the reactions of H 2 0 with three stabilized Criegee intermediates 
(stabCI-OO, stabCI-CH3-00 and stabChc-OO) from ozonolysis of limonene are shown in Schemes 
respectively. The structures of the intermediate adduct and transition states are denoted as M and TS, 
respectively. The HIS, H2_S and H3_S are used to distinguish the reactions of H 2 0 with stabCI-OO, 
stabCI-CH3-00 and stabChc-OO. The reactions of H2O with all the three stabilized Criegee 
intermediates are three-step reactions. For each reaction, a hydrogen-bond complex is initially formed. 
The further evolution of the hydrogen-bond complex formed occurs via two different 
reaction pathways: 

(1) Formation of OH radicals with water-catalyzed H migration. The water molecule assists the H 
migration toward the terminal oxygen of the COO moiety, which could lead to the formation of 
OH radicals. 

(2) Formation of a-hydroxy hydroperoxide. Cycloaddition of H2O to stabilized Criegee intermediate, 
where the oxygen of water is linked to the carbon atom of the stabilized Criegee intermediate 
while a hydrogen atom of water is transferred to the terminal oxygen of the COO unit. 

The reaction pathway (2) have different fates depending on the nature of the R and R f substitutes in 
subsequent reactions. For reaction between stabCl-00 and H 2 0, the M1H1_S (hydrogen-bond 
complex) is getting transformed via the two main reaction pathways (pathway (1) and pathway (2)) 
through the corresponding transition states (TSM11H1_S and TSM12H1_S), producing M2H1_S and 
M3H1_S, respectively. In pathway (1), M2H1S evolves via losing water molecule and the transition 
state (TSM4H1S), to produce OH and R1H1S radicals (pathway (aHl_S)). M3H1S in pathway (2) 
may have different degradation paths. M3H1_S can be transformed via the transition state 
(TSM31H1S) into OH and R2H1S radicals (pathway (blHl S)), limononic acid and H 2 0 
(pathway (b2Hl_S-b5Hl_S)) and limononaldehyde and H2O2 (via the transition states 
(TSM32H1_S-TSM35H1_S) (pathway (b6Hl_S)). Moreover, the pathway (b4Hl_S) was also found 
for the reaction between H 2 0 and stabCl-00 in the recent study of Wang [42]. 
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Scheme 2. Mechanistic diagram for the reaction between H2O and stabCI-00 arising from 
the limonene ozonolysis. 
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Scheme 3. Mechanistic diagram for the reaction between H 2 0 and stabCI-CH 3 -00 arising 
from the limonene ozonolysis. 




Limononaldehyde M5H2_S 



Int. J. Mol Sci. 2013, 14 



5788 



Scheme 4. Mechanistic diagram for the reaction between H2O and stabCbc-OO arising 
from the limonene ozonolysis. 
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In the reaction between stabCI-CH3-00 and H 2 0, the hydrogen-bonded complex M1H2_S evolves 
via their corresponding transition states (TSM11H2_S and TSM12H2S), leading to the formation of 
M2H2 S and M3H2S, respectively. In pathway (1), M2H2 S evolves via the detachment of the water 
molecule and the transition state (TSM4H2S), to produce OH and R1H2S radicals (pathway 
(aH2_S)). M3H2_S may follow three different reaction pathways in pathway (2). M3H2_S can be 
transformed into OH and R2H2_S radicals (pathway (blH2_S)) or via the transition states 
(TSM31H2 S and TSM32H2S), into H 2 0, R3H2S, R4H2S, limononaldehyde and H 2 0 2 . The 
reaction between stabCbc-OO and H 2 0 follows the same pathways as those for the reaction of 
stabCI-CH 3 -00 and H 2 0. The M1H3 S can evolve via their corresponding transition states 
(TSM11H3 S and TSM12H3 S) to formation of OH radical, R1H3S radical, R2H3 S radical, 
R3H3 S, R4H3 S, H 2 0, Keto-limonene and H 2 0 2 . 

The formations of limononic acid, limononaldehyde and keto-limonene are consistent with the 
experimental results for the reaction between three stabilized Criegee intermediates and H 2 0 in 
Limonene-ozone reaction in UNC smog chamber [32,43]. 

Although the important role of the reaction between Criegee intermediates from ozonolysis of 
limonene and H 2 0 for the formation of SO A and organic acids affecting both the global climate 
changes and public health is well established, the mechanism of these reactions remains poorly 
understood. Since the stability of the Criegee intermediates is moderately weak, the computational 
quantum chemistry is the most efficient for determining the geometric and electronic structures of 
these chemically activated complexes. Since reaction between the simple stabCH 2 00 and H 2 0 has 
been investigated in some detail [6,11,44-51] in the past, the mechanism of this reaction is left out of 
the scope of the present paper. 
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In the present paper, the reaction of H2O with three more complex stabilized Criegee intermediates 
(stabCI-OO, stabCI-CH 3 -00 and stabCbc-OO) from ozonolysis of limonene has been investigated in 
order to gain new insights of the oxidation mechanisms under atmospheric conditions. Ab initio and 
DFT methods have been employed to obtain geometries and energies of the transition states and 
subsequent degradation products. Reaction and activation energies for the reaction between the three 
stabilized Criegee intermediates and H 2 0 have been computed at different levels of theory, including 
the higher level CCSD(T)/6-31G(d) + CF. A thorough thermochemical analysis has been carried out 
and its results and implications for the ozonolysis of limonene have been discussed. 

2. Computational Details 

The computations were performed using GAUSSIAN 03 suite of programs [52] on the SGI ALTIX 
4700 supercomputer. The geometry optimization of all reactants, stabilized Criegee intermediates, 
transition states, and products was executed using Becke's three-parameter hybrid method employing 
the LYP correction function (B3LYP) with the split valence polarized basis set 6-31G(d,p) [10,53]. 
The stationary points were classified as minima in the case, when no imaginary frequencies were 
found, and as a transition state in the case, when only one imaginary frequency was obtained. In order 
to verify the transition states connecting reactants and products, the intrinsic reaction coordinate (IRC) 
analysis [54] at the B3LYP/6-31G(d,p) level of theory has also been applied to each transition state of 
every reaction. We just freely optimized the original three CI structure to get the minima conformer, 
then keeping the conformer for the succedent reactions, (see Figures 1-3). There exists difference 
between our optimized conformation and those in other studies. For example, in our study, we have 
chosen opposite orientation on stabCI-00 and stabChc-00 from what Baptista et al. did [55]. In 
general, we think our method of computational process is applicable for explanation of unanimous 
reaction between the same conformer of CI with water. The DFT structures were then used in the 
single-point energy calculations using frozen core second-order Nfeller-Plesset perturbation theory 
(MP2) and coupled-cluster theory with single and double excitations including perturbative corrections 
for the triple excitations (CCSD(T) [56] with various basis sets. The basis set effects on calculated 
energies for the reactions of stabilized Criegee intermediates with H 2 0 were corrected at MP2 level 
according to the developed method, which has been successfully applied for studying the complex 
reaction mechanisms and pathways of volatile organic compounds (VOCs) in the atmosphere [57]. A 
correction factor (CF) has been determined from the energy difference between the MP2/6-31G(d) and 
MP2/6-311++G(d,p) levels. Energies calculated at the CCSD(T)/6-31G(d) level of theory have been 
corrected using the aforementioned MP2 level corrections. The application of the CF has been 
validated in several studies of isoprene and limonene reactions initiated by NO3, OH and O3 [57-61]. 
For example, new important results have been obtained using the same method for the reaction of 
stabilized Criegee intermediates from the ozonolysis of limonene with sulfur dioxide in the recent 
paper [41]. These considerations lead us to conclude that the CCSD(T)/6-31G(d) + CF selected as a 
primary method and calculation scheme are appropriate for studying the reactions of VOCs with 
ozone or OH radical and description of the Criegee intermediates in the ozonolysis of VOCs 
reaction [35,37,47,62-66]. 
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Figure 1. Geometries of the stationary points in the stabCI-00 + H2O reaction obtained at 
the B3LYP/6-31G(d,p) level of theory. Bond lengths and intermolecular distances are 
given in A. 
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3. Results and Discussion 

3.1. Reaction Mechanism 

The reactions of H 2 0 and limonene stabilized Criegee intermediates (stabCI-OO, stabCI-CH 3 -00 
and stabChc-OO) occur via the six main reaction pathways. Each transition state has only one 
imaginary harmonic vibrational frequency and can be classified as the first-order saddle point. The 
values of imaginary frequencies for TSM11H1S, TSM12H1S, TSM11H2S, TSM12H2S, 
TSM11H3_S and TSM12H3_S transition states are 1516.76i, 434.95i, 1584.16i, 590.54i, 1586.3H and 
593.62i, respectively. 

Figure 1 presents the optimized geometries of the stationary points for the reaction of stabCI-00 
and H2O obtained at the B3LYP/6-31G(d,p) level of theory and the most important corresponding 
geometrical properties such as bond lengths and bond angles. Scheme 2 and Figure 1 show that the 
reactions of stabCI-00 with H 2 0 initially lead to the formation of the hydrogen-bond complex 
M1H1S. Then the M1H1S evolves via the transition state TSM11H1S (formation of OH radicals 
with water-catalyzed H migration path or reaction (1)), to the hydrogen-bond complex M2H1S, 
which can evolve into hydroperoxide M4H1_S and water. Finally, the cleavage of the O-OH bond in 
M4H1_S may lead to OH and R1H1_S radicals via the transition state TSM4H1_S along the reaction 
path (aHl_s). 

As seen from Figure 1, there are some major changes in main bonds. In particular, the 
hydrogen-bond M1H1_S complex includes a seven-membered ring. The distance between 04 and H9 
connecting the terminal oxygen atom (04) of the CI 0304 group belonging to stabCI-00 to the 
hydrogen atom (H9), is 1.837 A, while the length of the 08-H9 bond is 0.982 A. The length of H5-08 
bond, which connects the oxygen (08) of original water to the hydrogen atom (H5), is 2.363 A. The 
corresponding transition (TSM11H1_S) shows that the transfer of the hydrogen atom (H5) linked to 
the carbon atom (C2) to the water molecule can take place, and the transfer of one hydrogen atom (H9) 
from water to the terminal oxygen of the CI 0304 group can also occur at the same time. This 
indicates that the water molecule acts as a catalyst of the hydrogen migration. A comparison of the 
corresponding transition (TSM11H1_S) structure with M1H1_S shows that the C2-H5 bond length 
increases by 0.232 A to 1.331 A, while the H5-08 distance decreases by 1.031 A to 1.332 A. The 
08-H9 bond length increases by 0.166 A to 1.148 A, while the 04-H9 bond length decreases by 0.536 
A to 1.301 A. C2-H5 and 08-H9 distances in TSM11H1_S continue to increase, leading to both 
C2-H5 and 08-H9 bonds broken and the formation of 08-H5 and 04-H9 bonds. At the same time, 
the hydrogen-bond complex M2H1S is formed. The lengths of C2-H5, H5-08, 08-H9 and 04-H9 
bonds in M2H1_S are 2.567 A, 0.971 A, 1.804 A, 0.990 A and 1.446 A, respectively. A hydroperoxide 
M4H1_S and water are formed from the M2H1_S. The 03-04 and 04-H9 bonds in M4H1_S 
experiencing minor changes during this process are of 1.448A and 0.973 A, respectively. Finally, the 
O-OH bond in M4H1S is broken to form of 04H9 and R1H1S radicals via the TSM4H1S 
transition state. The 03-04 bond in the peroxide (2.049 A) in TSM4H1_S complex gets elongated by 
0.601 A compared to that in the M4H1_S. The reaction (1) may be considered as a possible significant 
source of atmospheric OH radicals. However, this pathway is a poor OH source. Because the vinyl 
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hydroperoxides (M4H1S) formed in this pathway is easily to collisionally stabilized under 
atmospheric condition and thus impeded OH formation [67,68]. 

M1H1_S can also evolve via the transition state TSM12H1_S (formation of a-hydroxy 
hydroperoxide or reaction (2)), into a a-hydroxy hydroperoxide M3H1_S. In this case, the formation of 
transition state TSM12H1_S is associated with the evolution of five-membered ring from the 
M1H1_S, where O8-H10 group and H9 atom belonging to the water move towards the CI atom and 
the terminal oxygen atom (04) of the CI 0304 group belonging to stabCI-OO, respectively. There are 
also significant changes in the length of the main bonds. In particular, the 04-H9 bond length (1.408 
A) in the corresponding transition (TSM12H1_S) complex decreases by 0.429 A and 08-H9 bond 
(1.103 A) increases by 0.121 A compared with those in M1H1_S. The distance between CI and 08 in 
the TSM12H1_S is 2.037 A. Then, the a-hydroxy hydroperoxide M1H1_S can be formed according to 
the following scheme: both CI -08 and 04-H9 distances in the TSM12H1_S continue to decrease, 
leading to 08-H9 bonds broken and the formation of both Cl-08 (1.402A) and 04-H9 (0.971A) 
bonds in the M3H1_S. 

The M3H1_S has six different reaction pathways (blHl_S-b6Hl_S). In the reaction pathway 
blHl S, M3Hl_s initially evolves via the transition state TSM31H1S. Then, the 03-04 bond 
(2.452 A) of TSM31H1_S is broken and the formation of R2H1_S and 04H9 radicals occurs. The 
blHl_S of reaction (2) may be considered as another possible source of atmospheric OH radicals. 

In the reaction pathway b2Hl_S/b3Hl_S, M3H1_S evolves via the transition state TSM32H1_S/ 
TSM33H1_S, yielding M5H1_S/M6H1_S. This process occurs according to the following scheme: the 
H5/H6 connected to C2 moves towards 04-H9 group, while H7 linked to CI migrates to C2 instead of 
H5/H6. After both 04-H5/04-H6 (0.967 A/0.965 A) and C2-H7 (1.093 A/1.906 A) bonds are formed, 
the formation of M5H1S/M6H1S is completed. The 03-04 bond increases to 2.274 A/2.268 A in 
TSM32H1S/TSM33H1S, and elongates to 3.148 A/2.923 A in M5H1S/M6H1S. Finally, the 
03-04 bond in the M5H1_S/M6H1_S is broken, leading to the formation of limononic acid 
A/limononic acid B and water. The reaction pathway b4Hl_S is similar to b2Hl_S. H7 connected to 
CI migrates to the 04-H9 group, forming 04-H7 bond. The M3H1_S transforms via the transition 
state TSM34H1S, into M7H1S. Finally, the limononic acid C and water are formed from M7H1S. 
The reaction pathway b5Hl_S is quite similar to b2Hl_S. However, in contrast to the b2Hl_S, H7 
connected to CI moves towards 03, while H10 bonded to 08 migrates to 04-H9 group. After 
O8-H10 bond is broken and both 03-H7 (0.989 A) and O4-H10 (0.975 A) bonds are formed, the 
formation M8H1_S is completed. Then, the 03-04 bond in the M8H1_S gets broken, leading to the 
formation of limononic acid D and water. In the reaction pathway b6Hl_S, M3H1_S evolves via the 
transition state TSM36H1_S, into M9H1S. This process occurs according to the following scheme: 
H10 connected to 08 moves towards 03 and the distance of O8-H10 increases by 0.720 A reaching 
1.687 A, while the Cl-03 bond is elongated by 1.980 A reaching 3.411 A in TSM36H1_S. 
TSM36H1_S then transforms into M9H1_S as the H10 connected to 08 continue to move towards 03 
and the distance of O8-H10 gets increased to 1.839 A, while the Cl-03 distance increases to 3.477 A. 
After the O3-H10 bond (0.983 A) is formed, the formation of M9H1_S is completed. Finally, the 
limono aldehyde and H2O2 are formed from M9H1S. The b6Hl_S can also be considered as a 
possible source of atmospheric H2O2. 
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Scheme 3 and Figure 2 show that the reaction of stabCI-CH3-00 with H2O exhibits the same 
behavior as stabCI-00 + H 2 0 reaction. At the first step, the reacting stabCI-CH 3 -00 and H 2 0 evolve 
into a hydrogen-bond M1H2_S complex of seven-membered ring. The M1H2_S finally transforms via 
the transition state TSM11H2_S (reaction (1)) into OH and R1H2_S radicals (pathway (aH2_S)) and 
via the TSM12H2 S (reaction (2)) into OH radical, R2H2S radical, H 2 0, R3H2S, R4H2S, 
limononaldehyde and H 2 0 2 (pathway (blH2_S-b3H2_S)). 



Figure 2. Geometries of the stationary points in the stabCI-CH3-00 + H 2 0 reaction 
obtained at the B3LYP/6-31G(d,p) level of theory. Bond lengths and intermolecular 
distances are given in A. 




Scheme 4 and Figure 3 also show that the reaction between stabCbc-OO and H 2 0 exhibits the same 
pattern and occurs according to the same scheme as those for the reaction between stabCI-CH 3 -00 and 
H 2 0. The reaction between stabChc-OO and H 2 0 leads to the formation of a hydrogen-bonded 
M1H3_S complex with seven-membered ring. Finally, the M1H3 S finally evolves via the transition 
states TSM11H3 S (reaction (1)) and TSM12H3 S (reaction (2)) into OH radical, R1H3S radical, 
R2H3 S radical, R3H3 S, R4H3 S, H 2 0, Keto-limonene and H 2 0 2 . 
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Figure 3. Geometries of the stationary points in the stabCIx-00 + H2O reaction obtained 
at the B3LYP/6-31G(d,p) level of theory. Bond lengths and intermolecular distances are 
given in A. 
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020103=125.1 
010304=111.7 
O3O4H7=100.3 
C2C1O3O4=-3.0 
J C10304H7=98. 




H5C1 02=1 07.1 
020103=119.3 
010304=117.8 
O3O4H7=100.0 
H506H7=73.6 
H5C201O3=-62.2 
C2C1 0304-1 .3 
C10304H7=72.3 
04H706H5=-5.8 



020103=116.3 
010304=87.2 
0304H7=134.1 
02010304=54.0 
C10304H7=164.9 




H5C1C2=98.5 
020103=122.1 
010304=116.0 
O3O4H7=100.0 
H506H7=86.6 
H5C2C103=-57.8 
C2C10304—1.2 A 
C10304H7=59.3^ 
04H706H5=-8.1 





020103=119.7 



.029 




H5C1 02=82.7 
020103=125.2 
q 74 010304=112.4 
.0304H7=99.3 
H5O6H7=90.5 
H5C201O3=-66.6 
02010304=0.6 
C10304H7=82.9 
04H706H5=-7.4 



C1O6H8=110.8 
010304=109.3 
0304H7=85.2 
H7O6H8=106.2 
H806C1 03=1 35.3 
09 6 7 06C10304=-51.3 
1112 01O3O4H7=48.1 
04H706H8=-114.4 

1.383 



010304=107.2 
O3O4H7=99.0 
C1O6H8=105.8 
030106=110.2 
090103=105.9 
C1C9H12=103.3 
01009011=109.9 
C10304H7=1 
03C106H8=67. 
09010306= 
H1 2C9C103: 



030106=114.2 



j TSM31H3_S 
154-2T ™1 fifll 




C106H8=86.3 
06H803=125.4 
■ 71 H80304=121.5 
0 971^O3O4H7=99.8 
A C1O6H8O3=-23.0 
^ O6H8O3O4=130.2 
H80304H7=161.1 




C106H8=117.6 
06H803=167.8 
H8O3O4=100.4 
O3O4H7=100.0 
C106H803=-5.3 
O6H8O3O4=102.3 
H80304H7=111.4 



Keto-limonene 



010304=92.9 
C1O6H8=94.0 
030106=113.8 
090103=87.2 
01 C9H1 2=89.0 
C10C9C11=114 
O3C1O6H8=30. 
C9C10306=-11 
H12C9C103=-1' 




01009011=119.3 
C10C9H12=120.2 
C10C9C11H12=16f 




3.2. Thermochemical Analysis 

Although the higher level CCSD(T)/6-31G(d) + CF method has been chosen as a principal 
computational method for the present study, in order to ensure the quality of the obtained results and to 
validate the obtained conclusions, reaction and activation energies for the reaction between stabCl-00 
and H2O with the Zero-Point Correction (ZPE) included were also calculated using several other 
methods used to study the similar systems in the past. The results of the present study are given in 
Table 1. As it may be seen from Table 1, the reaction and activation energies calculated at different 
levels are in qualitative agreement in most of the cases studied (not exceed 10 kcal mol 1 ). Only in a 
few cases, such as the energy of TSM31H1_S relative to M3H1_S, the values predicted 
B3LYP/6-31G(d,p), MP2/6-31G(d), MP2/6-311++G(d,p) and CCSD(T)/6-31G(d) are quite different 
from those given by the CCSD(T)/6-31G(d) + CF with the differences of 23.75 kcal mol" 1 , 
21.22 kcal mol 1 , 27.14 kcal mol 1 and 2.53 kcal mol 1 respectively, presumably most accurate method 
used in the present study. 
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Table 1. Reaction and activation energies (E) with Zero-Point Correction (ZPE) included 
(kcal mol 1 )) for the reaction between stabCl-00 and H2O a at different levels of theory. 



Compound 


Relative to 


^(B3LYP/ 


E(MP2/ 


E(MV2I 


^(CCSD(T)/ 


^(CCSD(T)/ 


6-31G(d,p)) 


6-31G(d)) 


6-311++G(d,p)) 


6-31G(d)) 


6-31G(d) + CF) 


stabCl-00 + H 2 0 




0.00 


0.00 


0.00 


0.00 


0.00 


A /T1 TT1 O 


stabCl-00 + 
H 2 0 


0 

-8.65 


-1 21 


-7.54 


-1.54 


-6.16 


TSM11H1S 


M1H1S 


13.27 


11.44 


8.52 


18.51 


16.68 


M2H1S 


TSM11H1S 


-34.17 


-34.57 


-26.94 


-35.09 


-35.50 


M4H1S + H 2 0 


M2H1S 


8.80 


5.64 


8.50 


7.91 


4.74 


TSM4H1S 


M4H1S 


24.29 


18.38 


29.60 


20.51 


14.60 


R1H1S+OH 


TSM4H1S 


-10.64 


-8.29 


-17.17 


-9.39 


-7.04 


TSM12H1S 


M1H1S 


14.96 


17.01 


13.62 


16.78 


18.83 


M3H1S 


TSM12H1S 


-50.15 


-50.35 


-41.65 


-47.92 


-48.12 


TSM31H1S 


M3H1S 


56.79 


54.26 


60.18 


35.57 


33.04 


R2H1S+OH 


TSM31H1S 


-12.35 


-12.64 


-22.88 


2.06 


1.78 


TSM32H1S 


M3H1S 


40.23 


32.82 


46.83 


49.85 


42.44 


M5H1S 


TSM32H1S 


-121.60 


-120.87 


-124.12 


-126.25 


-125.52 


Limononic acid 
A + H 2 0 


M5H1S 


7.41 


5.01 


6.31 


7.12 


4.71 


TSM33H1 S 


M3H1 S 


44.84 


35.62 


51.29 


54.09 


44.87 


M6H1 S 


TSM33H1 S 


-121.88 


-120.57 


-126.09 


-126.76 


-125.44 


Limononic acid B 
+ H 2 0 


M6H1S 


5.84 


3.56 


5.60 


5.53 


3.24 


TSM34H1 S 


M3H1 S 


53.34 


48.18 


51.67 


52.07 


46.91 


M7H1S 


TSM34H1S 


-130.73 


-131.72 


-126.97 


-124.94 


-125.93 


Limononic acid C 
+ H 2 0 


M7H1S 


8.06 


4.72 


7.79 


7.56 


4.22 


TSM35H1S 


M3H1S 


44.65 


42.72 


34.16 


43.70 


41.78 


M8H1S 


TSM35H1S 


-125.34 


-129.75 


-111.65 


-119.40 


-123.81 


Limononic acid 
D + H 2 0 


M8H1S 


8.44 


6.46 


7.78 


8.03 


6.05 


TSM36H1S 


M3H1S 


54.94 


55.16 


49.52 


52.33 


52.55 


M9H1S 


TSM36H1S 


-51.97 


-51.66 


-47.33 


-49.38 


-49.07 


Limononaldehyd 

e + H 2 0 2 


M9H1S 


8.58 


6.79 


8.12 


8.26 


6.47 



Optimized geometries, vibrational frequencies and ZPE obtained at the B3LYP/6-31G(d,p) level. 



As seen from Table 1 and Figure 4, the most accurate and higher level CCSD(T)/6-31G(d) + CF 
method show that the initial hydrogen-bond complex M1H1_S is 6.16 kcal mol 1 more stable than the 
separate stabCl-00 and H2O. For the reaction of stabCI-00 with H2O, the formation of a-hydroxy 
hydroperoxide (reaction (2)) (-35.45 kcal mol 1 ) is more favorable than the formation of OH radicals 
with water-catalyzed H migration path (reaction (1)) (-24.98 kcal mol 1 ) (see Figure 4) The energies 
for the reaction pathways aHl_S and blHl_S associated with the formation of OH radical are 
12.67 kcal mol 1 and 0.63 kcal mol 1 , respectively. Among the seven reaction pathways 
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(aHl_S-b6Hl_S), the pathways (aHl_S and blHl_S) for the formation of OH radical are by 
12.67 kcal mol 1 and 0.63 kcal moF 1 more stable than the separate stabCI-00 and H 2 0, the pathway 
(b6Hl_S) for formation of limonoaldehyde and H2O2 is by 25.49 kcal mol -1 more stable than the 
separate stabCI-00 and H2O, and the four pathways (b2Hl_S-b5Hl_S) for formation of limononic 
acid and water are by 1 10.25-1 13.82 kcal mol -1 more stable than the separate stabCI-00 and H 2 0. The 
differences in the reaction energies are close and agree within 3.57 kcal mol 1 for all the four pathways 
(b2Hl_S-b5Hl_S). 

Figure 4. stabCI-00 + H 2 0 reaction coordinates: relative energies of the stationary points 
located on the separate stabCI-00 and H 2 0 ground-state potential energy surface. The 
energy values are given in kcal mol -1 and were calculated using CCSD(T)/6-31G(d) + 
CF//B3LYP/6-31G(d,p) level. 



20 • 



-20 



-100 



-120 



12.6" 
10.53 > 



TSM12H1 S 



TSM36H1 S . 



sabCl-00 
+H,0 



M1H1_S/TSM11H1 s 



-6.16 



-24.98 
-35.45 



TSM34H1_S .. 

11.46 ;; # 
TSM33H1_S .*.V 
9.43 .7/ 



, M2H1JS "» 



TSM32H1 S 



6.99 



-TSM35H1 S 



6.33 



TSM31H1 S 



-2.40 R2Hl_S-OH 



\ V ' TSM4H1_S(H O) 
-5.63 



H O 



;M4H1_S(H 2 0){\ M9H1_S 

-20.23 i ' 

i -31.96 



-114.4" M~H1_S 
-116.01 M6H1_S 



-0.63 
R1H1_SOH(H,0) 



-12.67 



-25.49 



\\ M8H1_S . 117 .4 8 
M5H1_S -118.53 



. Limononaldehvde-H,0, 



L imononic acid OH,0 - 1 1 0. 25 

Limononic acid D~H,0 -111.43 
— — Limononic acid B-flO -11277 
Limononic acid A -H O -113.82 



With respect to the corresponding M1H1S, the activation energy for the formation of OH radicals 



with water-catalyzed H migration path (reaction (1)) (16.69 kcal moF 1 ) is smaller than that for the 
formation of a-hydroxy hydroperoxide (reaction (2)) (18.83 kcal mol 1 ), indicating that the formation 
of OH radicals with water-catalyzed H migration path (reaction (1)) is slightly more favorable than the 
formation of a-hydroxy hydroperoxide (reaction (2)). Thus, with the reaction between our optimized 
stabCI-00 (syn stabCI-OO) and water, this conclusion is opposite to that obtained in the earlier 
theoretical study of the reaction between isoprene stabilized Criegee intermediates and H2O by 
Anglada [37]. Since in Anglada's study [37], the energies for the reaction of water with all the eight 
isoprene stabilized Criegee intermediates have been calculated using G2(G2M-RCC5) method [38] 
with all stationary points optimized using the B3LYP/6-311+G(2d,2p) and neither comprehensive 
benchmarks for aforementioned method [38] nor a comparison between the prediction of the 
CCSD(T)/6-31G(d) + CF//B3LYP/6-31G(d,p) method and G2M-RCC5//B3LYP/6-311+G(2d,2p) 
methods are available at the present time, the activation energies between reaction (1) and reaction (2) 
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for stabCI-00 and H 2 0 obtained at the B3LYP/6-3 1 G(d,p)//B3LYP/6-3 1 G(d,p), 
MP2/6-3 lG(d)//B3LYP/6-3 lG(d,p), MP2/6-3 1 l++G(d,p)//B3LYP/6-3 lG(d,p) and CCSD(T)/6- 
31G(d)//B3LYP/6-31G(d,p) levels of theory in the present study were used to validate our conclusion. 
The comparison, which indicates that the reaction (1) is more favorable than reaction (2) for 
stabCI-00 and H 2 0 in all the cases studied (expect for the CCSD(T)/6-31G(d)//B3LYP/ 
6-31G(d,p)) method) largely confirms our conclusion. 

The activation energies of the six subsequent reaction pathways for the formation of a-hydroxy 
hydroperoxide (reaction (2)) are in the range of 33.04-52.55 kcal mol -1 with respect to corresponding 
M3H1_S, with the blHl_S as the most favorable pathway of the activation energy of 33.04 kcal mol 1 . 
The activation energies for the four pathways (b2Hl_S-b5Hl_S) are in the range of 
41.78-46.91 kcal mol 1 , and the differences between them are quite small (<3.578 kcal mol 1 ). 
Figure 4 illustrates the relative energies of the stationary points located on the singlet ground-state the 
separate stabCI-00 and H2O potential energy surface at the CCSD(T)/6-31G(d) + CF level of theory. 

Reaction and activation energies for the reaction between stabCI-CH 3 -00 and H 2 0 with the 
Zero-Point Correction (ZPE) included were calculated at different levels of theory. The results of the 
calculations are presented in Table 2. As it may be seen from Table 2 and Figure 5, the higher level 
CCSD(T)/6-31G(d) + CF predicts that the initial hydrogen-bond complex M1H2 S is 6.35 kcal mol" 1 
more stable than the separate stabCI-CH 3 -00 and H 2 0. For the reaction of stabCI-CH 3 -00 with H 2 0, 
the formation of a-hydroxy hydroperoxide (reaction (2)) (-32.17 kcal mol 1 ) is more stable than the 
formation of OH radicals with the water-catalyzed H migration path (reaction (1)) (-16.70 kcal mol 1 ). 
(See Figure 5) In the case of the four following reaction pathways (aH2_S-b3H2_S), the reaction 
energies with the respect to corresponding to the separated stabCI-CH 3 -00 and H 2 0, are 
-3.03 kcal mol 1 , -3.46 kcal mol 1 , -10.25 kcal mol 1 and -21.30 kcal mol 1 , respectively. The 
activation energies for the formation of OH radicals from M1H2 S via water-catalyzed H migration 
path (reaction (1)) (21.98 kcal mol 1 ) is larger than those for the formation of a-hydroxy 
hydroperoxide (reaction (2)) (15.53 kcal mol 1 ). This indicates that the formation of a-hydroxy 
hydroperoxide (reaction (2)) is more favorable than the formation of OH radicals with water-catalyzed 
H migration path (reaction (1)). The differences in the activation energies between reaction (1) and 
reaction (2) for stabCI-CH 3 -00 and H 2 0 obtained at the B3LYP/6-31G(d,p)//B3LYP/6-31G(d,p), 
MP2/6-3 lG(d)//B3LYP/6-3 lG(d,p), MP2/6-3 1 l++G(d,p)//B3LYP/6-3 lG(d,p) and CCSD(T)/6- 
31G(d)//B3LYP/6-31G(d,p) levels of theory also indicate that the reaction (2) is more favorable than 
reaction (1) for stabCI-00 and H 2 0 . This conclusion is in agreement with the previous theoretical 
study of the reaction between isoprene stabilized Criegee intermediates and H 2 0 by Anglada [37]. 

The activation energies of the three subsequent reaction pathways for the formation of a-hydroxy 
hydroperoxide (reaction (2)) are in the range of 35.63^44.04 kcal mol 1 with respect to corresponding 
M3H2_S, with the blH2_S as the most favorable pathway with the activation energy of 
35.63 kcal mol 1 . Figure 5 illustrates the relative energies of the stationary points located on the singlet 
ground-state the separate stabCI-CH 3 -00 and H 2 0 potential energy surface at the CCSD(T)/6-31G(d) 
+ CF level of theory. 
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Table 2. Reaction and activation energies (E) with Zero-Point Correction (ZPE) included 
(kcal mol 1 ) for the reaction between stabCl-CH3-00 and H2O a at different levels of theory. 



Compound 


Relative to 


£(B3LYP/ 


EG(MV2I 


E(MV2/ 


^(CCSD(T) 


^(CCSD(T)/ 


6-31G(d,p)) 


6-31G(d)) 


6-311++G(d,p)) 


/6-31G(d)) 


6-31G(d) + CF) 


stabCl-CH 3 -00 + H 2 0 




0.00 


0.00 


0.00 


0.00 


0.00 


M1H2S 


stabCl-CH 3 - 
00 + H 2 0 


-11.15 


-7.94 


-10.45 


-9.57 


-6.35 


TSM11H2S 


M1H2S 


19.60 


16.93 


14.68 


24.64 


21.98 


M2H2S 


TSM11 


-30.19 


-31.48 


-22.41 


-31.04 


-32.32 


M4H2S + H 2 0 


M2H2S 


8.74 


5.66 


8.71 


7.86 


4.79 


TSM4H2S 


M4H2S 


29.25 


25.80 


41.29 


20.30 


16.86 


R1H2S+OH 


TSM4H2S 


-17.41 


-15.58 


-28.27 


-9.81 


-7.98 


TSM12H2S 


M1H2S 


11.72 


13.61 


11.25 


13.64 


15.53 


M3H2S 


TSM12H2S 


-43.36 


-43.58 


-34.17 


-41.13 


-41.35 


R2H2S + OH 


M3H2S 


46.53 


43.33 


38.15 


38.84 


35.63 


TSM31H2S 


M3H2S 


45.47 


43.40 


37.90 


46.11 


44.04 


R3H2S + R4H2S + 
H 2 0 


TSM31H2S 


-20.33 


-22.99 


-21.49 


-19.45 


-22.11 


TSM32H2S 


M3H2S 


44.87 


43.12 


36.08 


44.51 


42.76 


M5H2S 


TSM32H2S 


-41.41 


-39.53 


-35.87 


-41.49 


-39.62 


Limononaldehyde + 
H 2 0 2 


M5H2S 


9.16 


8.17 


8.14 


8.72 


7.73 



Optimized geometries, vibrational frequencies and ZPE obtained at the B3LYP/6-31G(d,p) level. 



Reaction and activation energies for the reaction between stabCIx-00 and H2O with the Zero-Point 
Correction (ZPE) included calculated at different levels of theory are shown in Table 3. Table 3 and 
Figure 6 demonstrate that the values of reaction and activation for the reaction between stabChc-OO 
and H2O correlate with those of reaction and activation energies for the reaction of stabCI-CH3-00 
and H2O. The formation of a-hydroxy hydroperoxide (reaction (2)) is more favorable than the 
formation of OH radicals with water-catalyzed H migration path or reaction (1). The differences in the 
activation energies between reaction (1) and reaction (2) for stabChc-OO and H2O obtained at the 
B3LYP/6-3 lG(d,p)//B3LYP/6-3 lG(d,p), MP2/6-3 lG(d)//B3LYP/6-3 lG(d,p), MP2/6-3 1 1++ 
G(d,p)//B3LYP/6-3 1 G(d,p) and CCSD(T)/6-3 1 G(d) //B3LYP/6-3 1 G(d,p) levels of theory 
(4.06 kcal mol 1 , 0.12 kcal mol 1 , -0.73 kcal mol 1 and 7.16 kcal mol 1 , respectively) also indicate 
(with exception for MP2/6-311++G(d,p)//B3LYP/6-31G(d,p) method) that the reaction (2) is likely 
more favorable than reaction (1) for stabChc-OO and H 2 0. These findings are consistent with the 
previous theoretical study of the reaction between isoprene stabilized Criegee intermediates and H2O 
by Anglada [37]. Figure 6 illustrates the relative energies of the stationary points located on the singlet 
ground-state the separate stabChc-OO and H 2 0 potential energy surface at the CCSD(T)/6-31G(d) + 
CF level of theory. 
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Figure 5. StabCI-CH3-00 + H2O reaction coordinates: relative energies of the stationary 
points located on the separate stabCI-CH 3 -00 and H 2 0 ground-state potential energy 
surface. The energy values are given in kcal mol 1 and were calculated using CCSD(T)/6- 
31G(d) + CF//B3LYP/6-31G(d,p) level. 



20 
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© o-| 
E 
53 



u 

o 
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20- 



-30- 



-40 



TSM11H2 S . 



15.62 



TSM12H2 S • 918 



stabCl-CH 3 -00 
+HO 



-6.35 ♦ 
M1H2 S 



11.87 

TSM31H2S, 

TSM32H2S;,' 1059 
3.46 



R2H2 S+OHV 



-11.91 



M2H2_S H M4H2_S(H 2 0) 
'HO 



-16.70 



-32.17 : 

M3H2 S 



TSM4H2_S(H,0) 



4.94 



RlH2_S+OH(H O) 



-3.03 

\R3H2 S+R4H2 S+H O 



-10.25 



Limononaldehvde+H O 



-21.30 



-29.03 



M5H2 S 



Table 3. Reaction and activation energies (E) with Zero-Point Correction (ZPE) included 
(kcal mol -1 ) for the reaction between stabCbc-00 and H 2 0 a at different levels of theory. 



Compound 


Relative to 


^(B3LYP/ 


E(MF2/ 


E(MF2I 


^(CCSD(T) 


^(CCSD(T)/ 


6-31G(d,p)) 


6-31G(d)) 


6-311++G(d,p)) 


/6-31G(d)) 


6-31G(d) + CF) 


stabClx-00 + H 2 0 




0.00 


0.00 


0.00 


0.00 


0.00 


M1H3S 


stabClx-00 + 
H 2 0 


-10.79 


-7.67 


-10.22 


-9.32 


-6.19 


TSM11H3S 


M1H3S 


15.56 


12.93 


11.36 


20.55 


17.92 


M2H3S 


TSM11H3S 


-30.54 


-31.66 


-22.01 


-31.27 


-32.39 


M4H3S + H 2 0 


M2H3S 


9.31 


6.32 


8.79 


8.35 


5.36 


TSM4H3S 


M4H3S 


28.12 


24.23 


34.67 


24.14 


20.24 


R1H3S+OH 


TSM4H3S 


-12.80 


-11.39 


-18.84 


-10.79 


-9.38 


TSM12H3S 


M1H3S 


11.50 


12.81 


12.09 


13.39 


14.70 


M3H3S 


TSM12H3S 


-42.89 


-42.86 


-33.40 


-40.72 


-40.69 


R2H3S + OH 


M3H3S 


45.32 


42.42 


35.58 


37.28 


34.38 


TSM31H3S 


M3H3S 


41.90 


38.07 


34.61 


43.19 


39.36 


R3H3S + R4H3S + 
H 2 0 


TSM31H3S 


-19.76 


-20.40 


-24.55 


-20.32 


-20.95 


TSM32H3S 


M3H3S 


42.47 


41.24 


34.04 


42.27 


41.04 


M5H3S 


TSM32H3S 


-40.78 


-38.49 


-36.13 


-41.14 


-38.86 


Keto-limonene + 
H 2 0 2 


M5H3S 


8.20 


6.41 


7.99 


7.88 


6.10 



a Optimized geometries, vibrational frequencies and ZPE obtained at the B3LYP/6-31G(d,p) level. 
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As it may be seen from Tables 1-3, we can conclude that the energy of the initial hydrogen-bond 
complex (M1H1_S, M1H2_S, M1H3_S) are 6.16 kcal mol 1 , 6.35 kcal mol 1 and 6.19 kcal mol 1 more 
stable than the corresponding separate stabilized Criegee intermediates (stabCl-OO, stabCI-CH3-00, 
stabCLr-OO) and H 2 0 at the CCSD(T)/6-31G(d) + CF level of theory. For the formation of OH 
radicals via the water-catalyzed H migration path (reaction (1)), the activation energies with the respect 
to hydrogen-bond complex (M1H1_S, M1H2_S, M1H3_S) are 16.68 kcal mol" 1 , 21.98 kcal mol" 1 , and 
17.92 kcal mol" 1 , respectively, with the reaction (1) for the reaction of stabCl-00 and H2O as the most 
favorable pathway (16.68 kcal mol 1 ). For the formation of a-hydroxy hydroperoxide (reaction (2)), 
the activation energies are in the range of 14.70-18.83 kcal mol 1 with respect to corresponding 
hydrogen-bond complex (M1H1S, M1H2 S, M1H3 S), with the reaction (2) for stabCIx-00 + H 2 0 
reaction as the most favorable pathway (14.70 kcal mol" 1 ). 

Figure 6. stabCIx-00 + H 2 0 reaction coordinates: relative energies of the stationary points 
located on the separate stabCLx-00 and H 2 0 ground-state potential energy surface. The 
energy values are given in kcal mol 1 and are calculated using CCSD(T)/6-31G(d) + 
CF//B3LYP/6-31G(d,p) level. 
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4. Conclusions 



In the present study, several important aspects of the gas-phase reaction of the stabilized Criegee 
intermediates and H2O have been investigated. 

The present study leads us to the following conclusions: 

(a) The reaction between the stabilized Criegee intermediates (stabCI-OO, stabCI-CH 3 -00 and 
stabChc-OO) and H2O is the three-step reaction. At the first stage, the formation of a 
hydrogen-bonded complex occurs. At the second stage, the reaction can proceed via the 
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following two reaction pathways: (1) formation of OH radicals with water-catalyzed H 
migration; and (2) formation of a-hydroxy hydroperoxide. The formation of a-hydroxy 
hydroperoxide (reaction (2)), the reaction of stabCI-00 and H 2 0 occurs via six different 
degradation pathways, while the reaction of stabCI-CH3-00 and stabCbc-00 with H2O occurs 
via three different reaction pathways. 

(b) The formation of OH radicals with water-catalyzed H migration path (aHl_S, aH2_S and 
aH3_S) and formation of a-hydroxy hydroperoxide (reaction pathways blHl_S, blH2_S and 
blH3_S) may be considered as possible sources of OH radicals in the Earth's atmosphere. The 
formation of a-hydroxy hydroperoxide (reaction pathways b6Hl_S, b3H2_S and b3H3_S) can 
be considered as possible sources of atmospheric H2O2. 

(c) The CCSD(T)/6-31G(d) + CF activation energies are in the range of 14.70-21.98 kcal mol" 1 
with respect to corresponding hydrogen-bond complexes (M1H1_S, M1H2_S, M1H3_S) 
between stabilized Criegee intermediates and H2O. The formation of a-hydroxy hydroperoxide 
(reaction (2)) for the reaction of stabCbc-OO and H 2 0 is the most favorable pathway with the 
activation energy of 14.70 kcal mol 1 . 

(d) For the reaction of stabCI-00 and H2O, with the reaction between our optimized stabCI-00 
(syn stabCI-OO) and water, the formation of OH radicals with water-catalyzed H migration path 
or reaction (1) is more favorable than the formation of a-hydroxy hydroperoxide or reaction (2). 
This conclusion is opposite to that obtained in the previous theoretical study of the reaction 
between isoprene stabilized Criegee intermediates and H 2 0 by Anglada. For the reaction of 
stabCI-CH 3 -00 and stabCbc-OO with H 2 0, the formation of a-hydroxy hydroperoxide (reaction 
(2)) is more favorable than the formation of OH radicals with water-catalyzed H migration path 
or reaction (1). This conclusion is consistent with results of the previous theoretical study of the 
reaction between isoprene stabilized Criegee intermediates and H 2 0 by Anglada. 
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